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Abstract 

We study the neutralino dark matter phenomenology in the context of effective field 
theories derived from the weakly-coupled heterotic string. We consider in particular 
direct detection and indirect detection with neutrino telescopes rates. The two cases of 
moduli dominated and dilaton dominated SUSY breaking lead to completely different 
phenomenologies. Even if in both cases relic density constraints can be fulfilled, moduli 
domination generically leads to detection rates which are much below the present and 
future experimental sensitivities, whereas dilaton domination gives high detection rates 
accessible to the next generation of experiments. This could make dark matter searches 
an alternative way to constrain high energy fundamental parameters. 
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Conclusio 



1 Introduction 



There exists a large collection of measurements providing convincing evidence in favor 
of the existence of cold dark matter in the universe pQ. But the exact nature of this 
dark matter is still an open question. One of the most promising and best motivated 
candidates is a Weakly Interacting Massive Particle (WIMP). Direct detection via its 
rare scattering with a nucleus in a detector, or indirect detection via its annihilation 
after gravitational storage in a massive body provide two possible experimental strate- 
gies. 

In the framework of supersymmetry (SUSY), most extensions of the Standard 
Model (SM) predict a massive neutral weakly interacting particle in the form of a 
neutralino (xi = x)- Moreover, in the simplest versions of SUSY models such as the 
minimal supergravity model (mSUGRA),this particle is predicted to be stable. Recent 
works have constrained a large part of the parameter space available in mSUGRA [2] 
in constrained or unified versions of the model. For small values of Mq, we can even 
find in the literature some strong consequences on the limit of the neutralino mass : 
< 500 GeV [H] if we took into account all the recent accelerator analysis, but this can 
be strongly evaded when one allows high values of Mq and m^2 HI- Up to now, several 
works have generalized these simple models to a non-universal framework in the higgs 
sector El 13 IE] j the gaugino sector [HI IH1 HQ] or the sfermion sector . Such analyses 
show that neutralino dark matter searches are sensitive to the spectrum of supersym- 
metric particles, but the direct connection with the supersymmetry-breaking sector is 
not made because it appears only through the parameters of the effective theory, the 
so-called soft supersymmetry breaking parameters. 

If one wants to make explicit the supersymmetry breaking mechanism, one must 
identify the type of breaking as well as the nature of the mediation between the super- 
symmetry breaking sector and the observable sector. A standard example is gaugino 
condensation in a hidden sector which interacts only gravitationally with the sector 
of quarks and leptons. Since this involves gravitational interactions, it is natural to 
consider such a model in a string context and thus to include all fields associated with 
the gravitational sector, in particular the dilaton (whose vacuum expectation value 
determines the magnitude of the string coupling) and moduli fields (whose expecta- 
tion values determine the size of the compact manifold). Both types of fields play an 
important role in supersymmetry breaking. The most elaborated class of models and 
probably the most realistic one from a phenomenological point of view is associated 
with the weakly coupled heterotic string compactified on a Calabi-Yau manifold or an 
orbifold. 

Recently, the full one loop soft supersymmetry breaking terms in a large class of 
superstring effective theories have been calculated J2j based on orbifold compactifi- 
cations of the weakly-coupled heterotic string (including the so-called anomaly me- 
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diated contributions). Such models yield specific non-universalities which make their 
phenomenology significantly different from the minimal supergravity model. The pa- 
rameter space in this class of models has already been severely constrained by taking 
into account accelerator and relic density constraints |13l I14j or benchmark models at 
the Tevatron [151117)] . 

In what follows, we take this specific class of models to discuss how direct and 
indirect neutralino detection depend on the properties of supersymmetry breaking. 
Indeed, supersymmetry breaking induced by gaugino condensation may be transmitted 
through the auxiliary field vacuum expectation values of the compactification moduli 
or of the dilaton. We will see that the two corresponding cases have a very different 
behavior regarding dark matter detection. 

2 Theoretical framework 

2.1 Structure of heterotic orbifolds models at one loop 

The task of string phenomenology is to make contact between the high energy string 
theory, and the low energy world. For this purpose, we need to build a superstring 
theory in four dimensions, able to give us the Standard Model gauge group, three 
generations of squarks, and a coherent mechanism of SUSY breaking. We will set here 
our analysis in the framework of orbifold compactifications of the heterotic string within 
the context of the supergravity effective theory. We concentrate on those models where 
the action is dominated by one loop order contributions to soft breaking terms. The key 
property of such models is the non-universality of soft terms, consequence of the beta- 
function appearing in the superconformal anomalies. This non-universality gives a 
particular phenomenology in the gaugino and the scalar sector, modifying considerably 
the predictions coming from mSUGRA. In fact, these string-motivated models show 
new behavior that interpolates between the phenomenology of unified supergravity 
models (mSUGRA) and models dominated by the superconformal anomalies (AMSB). 
The constraints arising from accelerator searches and relic density have been already 
studied in jTH] . It is thus interesting, to see to which extent direct and indirect detection 
of dark matter will be able to bring us extra information on the models. 

We provide a phenomenological study within the context of orbifold compactifi- 
cations of the weakly-coupled heterotic string, where we distinguish two regimes. In 
the first one, the SUSY breaking is transmitted by the compactification moduli T a , 
whose vacuum expectation values determine the size of the compact manifold. Generic 
(0,2) orbifold models contain three T a moduli fields. We considered a situation in 
which only an " overall modulus T" field contributes to SUSY-breaking. The use of an 
overall modulus T is equivalent to the assumption that the three T a fields of generic 
orbifold models have similar contributions to SUSY-breaking. This is expected in the 
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absence of some dynamical effect that would strongly discriminate between the three 
moduli. In our second example, it is the dilaton field S present in any four-dimensional 
string (whose vacuum expectation value determine the magnitude of the unified cou- 
pling constant gsTR at the string scale), that transmits, via its auxiliary component, 
the SUSY breaking. We work in the context of models in which string nonperturba- 
tive corrections to the Kahler potential act to stabilize the dilaton in the presence of 
gaugino condensation |17l ITS] . The origin of the soft breaking terms are completely 
different in the two scenarii. Some are coming from the superconformal anomalies 
and are non-universal (proportional to the beta-function of the Standard Model gauge 
groups), others are generated in the hidden sector (from Green-Schwarz mechanism 
or gaugino condensation) and are thus universal. This mixture between universality 
and non-universality gives the richness of the phenomenology in this type of effective 
string models and confirms the interest of non-universal studies in the prospect of su- 
persymmetric dark matter detection, the non-universality being in this case connected 
with the basic properties of the model. 



2.2 The moduli dominated scenario 

In the moduli dominated scenario, the one loop order supersymmetric SUSY breaking 
terms at GUT scale can be written [T2 ] IT9| 120]: 
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(2.1) 
(2.2) 
(2.3) 



where M a and M, are the soft masses for the gauginos and scalars and Ai, the trilinear 
coupling. b a is the beta-function coefficient for the gauge group G a : 
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(2.4) 



where C a , C l a are the quadratic Casimir operators for the group G a in the adjoint 
representation and in the representation of the field i respectively. F s and F T are 
the auxiliary fields for the dilaton and the Kahler modulus, respectively, M is the 
supergravity auxiliary fields whose vacuum expectation value (vev) determines the 
gravitino mass 1713/2 = —\M, and 8qs is the Green-Schwarz coefficient which is a 
(negative) integer between and —90. The function G 2 {T, T) is proportional to the 
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Eisenstein function and vanishes when T is stabilized at one of its two self-dual points. 
From Ea, ()2.1|) . it follows that when the moduli are stabilized at a self dual point, 
only the second term contributes to gaugino masses. This is precisely the "anomaly 
mediated" contribution. The loop contributions have been computed using the Pauli- 
Villars (PV) regularization procedure. The PV regular fields mimic the heavy string 
modes that regulate the full string amplitude. The phenomenological parameter p 
which represents the effective modular weight of the PV fields is constrained to be 
no larger than 1, though it can be negative in value. Thus the scalar squared mass 
for all matter fields is in general non-zero and positive at one loop (only the Higgs 
can have a negative running squared mass). The limiting case of p = 1, where the 
scalar masses are zero at one loop level and for which we recover a sequestered sector 
limit, occurs when the regulating PV fields and the mass-generating PV fields have 
the same dependence on the Kahler moduli. Another reasonable possibility is that the 
PV masses are independent of the moduli, in which case we would have p = 0. and 7, 
is related to the anomalous dimension through y 1 - = ^i5j. (see |12H13lfT3| for notations 
and conventions) 

We clearly see in these formulae the competition between universal terms and non- 
universal ones. The scalar mass terms are all non-universal and proportional to their 
anomalous dimension 7, and thus loop suppressed. The Green-Schwarz mechanism 
generates universal breaking terms for the gauginos (proportional to 5qs) whereas 
superconformal anomalies introduce non— universal contributions (proportional to b a ^. 
The nature of the neutralino thus depends mainly on the value of the Green-Schwarz 
counterterm 5qs, whereas the mass scale is the gravitino mass 777.3/2- 

2.3 The dilaton dominated scenario 

We turn now to a scenario where the dilaton is the primary source of supersymmetry 
breaking in the observable sector. It is well known that if we use the standard Kahler 
potential derived from the tree level string theory, it is very difficult to stabilize the 
dilaton at acceptable weak-coupling values. We postulate in our study nonperturbative 
correction of stringy origin to the dilaton Kahler potential. In that case, one condensate 
can stabilize the dilaton at weak coupling while simultaneously ensuring vanishing 
expectation values at the minimum of the potential. The key feature of such models 
is the deviation of the dilaton Kahler metric from its tree level value. If we imagine 
the superpotential for the dilaton having the form W(S) oc e~ 3S ^ b+ , with b + being 
the largest beta-function coefficient among the condensing gauge groups of the hidden 
sector, then we are led to consider the phenomenology of models given by the following 
pattern of soft supersymmetry breaking terms |12 [ 1191 OH] : 
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+ h.c. 



where /i uv is an ultraviolet regularization scale (of the order of the string scale Mstr) 
and fj,R the renormalization scale (taken at the boundary value of Mqut) 1 - Moreover 



F s = v / 3m 3/2 ( J PC ss -)- 1 / 2 , K sS = d s dgK. 



(2- 
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with <7str the unified constant at -Mstr (see ^1 El Q5| f° r the notations and conven- 
tions) to ensure a vanishing vacuum energy in the dilaton-dominated limit. 

The phenomenology of a dilaton-dominated scenario is completely different from a 
moduli-dominated one. If we look at formulae (|2.5|) and (|2,7|) . it is clear that we are 
in a domain of heavy squarks and sleptons (of the order of magnitude of the gravitino 
mass) and relatively light gauginos. Indeed, the factor 6 + , as it contains a loop factor, 
can suppress the magnitude of the auxiliary field F s relative to that of the supergravity 
auxiliary field M through the relation (|2.8|) . The resulting gaugino soft breaking terms 
are less universal for low values of b + . 



1 For simplicity, we assume here that Mstr 
and appears in a loop factor). 



jU w ~ I^R ~ Mqxjt (the corresponding error is logarithmic 
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3 Supersymmetric dark matter phenomenology 



3.1 Relic density 

We recall that, in a general supersymmetric model, the neutralino mass matrix reads 
in the bino, wino, higgsino basis (B, W 3 , H®, H 2 ) : 



M 



N 



/ Mi — mz cos (3 sin W mz sin (3 sin 6 W \ 

M 2 mz cos P cos W —mz sin (3 cos 9 W 

-mz cos f3 sin 9 W mz cos /? cos 9 W — // 

y mz sin /? sin 0^ — sin /3 cos 9 W —\i / 

(3.10) 

Thus the lightest neutralino x is generically a superposition of these states : 

X = z X)1 B + z x . 2 W + z x ,3#i + z xA H 2 (3.11) 

In a large parameter space of the mSUGRA model, the x is mainly bino-like because 
of the renormalization group evolution of vr? H down to the electroweak scale. Indeed, 
let us recall the radiative electroweak symmetry breaking relation 

where &m 2 H ^ and &m 2 Hd represent the one loop tadpole corrections to the running Higgs 
masses m 2 H ^ and m 2 Hd |211 1221 123j . In (|3.12|) all the parameters are running and set at 
the minimization scale. We see that, as is becoming negative and relatively large 
in absolute value at low energy, which sets the higgsino mass scale, becomes large, 
in fact larger than the gaugino mass terms Mi and M 2 . Since in mSUGRA models, 
Mi = x )2 are linked at the weak scale by the well- known relation Mi = | tan 2 9\yM 2 ~ 
0.5M 2 (with some corrections at 2-loop order running), the x aligns along its bino 
component. When m 2 H (negative) decreases in absolute value, [i is then smaller, 
increasing the higgsino content of the neutralino. In the | m 2 H ^ \ decreasing direction, this 
leads successively to a mixed bino-higgsino neutralino and then an higgsino LSP and 
finally the no EWSB boundary when [i = 0. This happens typically in the hyperbolic 
branch (focus point) mSugra region thanks to the heavy scalar scale. 

Different processes lead to a cosmologically favoured neutralino relic density. For 
a bino neutralino one needs light sfermions, f(i) coannihilation or annihilation into 
pseudo-scalar A. If the lightest neutralino has a dominant wino component, the relic 
density drops down because of efficient annihilations into gauge bosons as well as strong 
XXi coannihilations. For a non negligible higgsino component neutralino annihilates 
into gauge bosons or ti and relic density is also decreased by xxt an d XX 2 coannihila- 
tions. 
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3.2 Direct detection 



Direct detection consists of measuring the energy deposited in a low background de- 
tector by the recoil of a nucleus from its elastic scattering with a Weakly Massive 
Interacting Particle (WIMP) In our case, the best WIMP candidate is the light- 
est neutralino x- The rates follow the neutralino-proton spin-independent elastic cross 
sections (<7*™,) or spin-dependent one (cr^!™), function on the target nucleus spin 

E3EE1- ' 

&x-p * s essen tially driven by first generation scalar quark (uj, di) exchanges or 

neutral Higgs (h, H) (xQ XQ) an d the spin-dependent one crf^™ by first generation 

squark and Z exchanges (xq — ^ XQ)- The processes involving Z-boson exchange being 
completely dependent on the neutralino higgsino fraction. 

The main impact of the null searches, particularly at LEP, is in the increase in the 
lower limit to the LSP mass, as well as the rest of the sparticle spectrum. 

In the mSUGRA case, the bino-like nature of the lightest neutralino x implies a 
highly suppressed scalar cross section via heavy neutral higgs exchange because of 
the low couplings x^X (proportional to the product of their higgsino and gaugino 
component) at moderate value of tan (3. If we look at the spin-dependent cross section, 
dominated by Z exchange, it becomes much larger when % is mostly higgsino because of 
the enhancement due to the x%X coupling (proportional to the square of the higgsino 
components). 

The real possibilities of an enhancement of u^-p are thus in the high tan j3 regime 
or in models that predict low first generation squarks masses for moderate gaugino 
masses, like some of the string inspired models that we study in this paper. Another 
possibility is to find regions of parameter space that increase the higgsino component of 
the lightest neutralino, enhancing in the same way its coupling. This happens for large 
Mq in mSUGRA along the "no EWSB" boundary. We will see that for a large class of 
heterotic orbifold models, all these constraints can be achieved. At the opposite, some 
of them will exhibit a complete depletion of the scalar cross-section, far below the 
sensitivity of the next generation of detectors. In this sense, direct detection of dark 
matter can become an important tool in the effort of discriminating string inspired 
models and constraining their parameter space. 

The current experimental status may be briefly summarized as follows. Although 
one of the current experiments, the DAMA collaboration 2x] claimed an evidence and 
gave a determination of the allowed maximum-likehood region in the WIMP-mass 
and WIMP-nucleon cross section of 1CP 6 — 10 _5 pb for a WIMP's mass between 30 
and 270 GeV, other experiments exclude almost (CDMS [H], EDELWEISS |2S|) or 
completely (ZEPLIN I [HJ) the DAMA region. But many new or upgraded versions 
of direct detection experiments will soon reach a significantly improved sensitivity for 
WIMP detection (EDELWEISS II [SO], ZEPLIN(s) [SI]). Our study will be placed 
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in the light of this next generation of detectors to see whether a SUSY dark matter 
candidate could be directly detected in the next years, and how this would constrain 
some of the fundamental SUSY breaking terms, and, as a consequence, the parameters 
of the more fundamental string theory. 

3.3 Neutrino indirect detection 

Neutralinos can also be gravitationally captured in massive astrophysical bodies like 
the Sun or the Earth by successive diffusions on their nuclei leading to a trapped neu- 
tralino population at their center. Then neutralinos can annihilate and the annihilation 
products, essentially gauge bosons and heavy quarks, decay emitting neutrinos. After 
conversion into muons through the Earth, these neutrinos can be observed by neutrino 
telescopes collecting Cerenkov light of induced muons traveling in water or in ice. The 
annihilation rate depends both on the capture rate (c*^, " depending on the target 
nucleus spin) and on the neutralino annihilation cross section |25| . The capture rate 
in the Earth depends on a^: p because of the zero spin of the iron nucleus. This leads 
to muon fluxes far beyond reach of detection [22]. In the case of the Sun, capture 
rates are enhanced by Z exchange in cr^ p thanks to the non zero spin of the hy- 
drogen. Fluxes are then maximized for a substantial higgsino fraction. Furthermore, 
mixed higgsino neutralino states annihilate into W + W~, ZZ or ti leading to more 
energetic neutrinos/muons than other annihilation channels. This happens along the 
"noEWSB" boundary, where (thanks to rrijj running) the neutralino gets a dominant 
higgsino fraction. This is the case in mSUGRA for high values of Mq |33l I32| and can 
be strongly favored with non universal gaugino masses by decreasing M^\gut |Hj- 

Some experiments have constrained these fluxes (Macro 34 , Baksan |35j . Super- 
Kamiokande [3*H]). but future neutrinos telescopes like Antares (3Z1 and Icecube j3H] will 
be much more efficient. They will improve current sensitivities of order 5.10 3 /i km -2 yr~ 
to 10 3 - 10 2 n km" 2 yr 1 on muon fluxes coming from neutralino annihilations in the 
centre of the Sun. We will compare our predictions for fluxes coming from the Sun 
with both current and future sensitivities. 

As we will see below, in the models that we are considering in this work, the muon 
fluxes coming from the Sun are very dependent on the nature of SUSY breaking, i.e. 
whether it is dominated by moduli or dilaton F terms. In this way, a neutralino dark 
matter signal in a neutrino telescope would provide key information on the nature of 
SUSY breaking and on the fundamental underlying theory. 

Works on prospects for direct and/or neutrino indirect detection of neutralino dark 
matter in mSUGRA/CMSSM or non universal frameworks can be found in |H 151 13H1 

El El El E01 EH E21 SSI G2H1 E21 [TT] . 
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3.4 Tools and experimental constraints 



In this section we describe the tools that we have used and the various constraints that 
we have imposed to obtain a correct phenomenology at the electroweak scale for our 
analysis. 

We have used the Fortran code SuSpect2 [44] to solve the renormalization group 
equations (RGEs) for the soft supersymmetry breaking parameters between the high 
energy boundary scale /i uv and the scale given by the Z-boson mass (electroweak scale) . 
While the initial scale u vv should itself be treated as a model-dependent parameter, 
for our purposes we have chosen for u vv the scale of grand unification M GUT . We use 
tan (3 and the sign of the supersymmetric [i parameter in the superpotential as free 
parameters, defined at the low-energy (electroweak) scale. 

The magnitude of the u parameter is determined by imposing electroweak symmetry 
breaking (EWSB) at the usual scale (mj m^J 1 / 2 (33 IS]- The one-loop corrected u is 
obtained from the condition (|3.12j) . 

The soft supersymmetry breaking parameters at the weak scale are then passed on 
to the C code micrOMEGAs [Uj to perform the calculation of physical masses for the 
superpartners and various indirect constraints, to be described below. 

We note that the value of the u parameter is fundamental in the analysis of astropar- 
ticle processes, because it determines the nature, mass and couplings of the lightest 
neutralino x which we require to be the LSP. We then estimated the detection rates 
using the DarkSusy package j3H] • 

The remaining parameter space is further reduced by limits on superpartner and 
Higgs masses from various collider experiments. We take the most recent bounds given 
by the different experiments of the LEP Working Group |49[ I50j . Concerning the light 
CP-even neutral Higgs mass (rrih), we assume that a 95% confidence level (CL) lower 
limit on is set at 111.5 GeV. The search for an invisibly decaying Higgs boson in hZ 
production has allowed a 95% CL lower limit on to be set at 114.5 GeV, assuming 
a production cross section equal to that in the Standard Model and a 100% branching 
fraction to invisible decays [HJ- We believe the value of 113.5 GeV will serve as a 
good mean. Concerning the chargino limit, we take 103.5 GeV, bearing in mind that 
in some degenerate cases and for light sleptons the limit can go down to 88 GeV |52j . 
For the squark sector the limit of 97 GeV jSH] is used. For all mass bounds we should 
keep in mind that experimental limits are always given in the context of a particular 
SUSY model framework which is not generally a string motivated one. The bounds we 
use could possibly be weakened in some cases. 

Various non-collider observations can be used to further reduce the allowed param- 
eter space of the loop-dominated orbifold models that we consider. We will focus our 
attention on the three sets of data that are the most constraining for these models: the 
density of relic neutralino LSPs, the branching ratio for decays involving the process 
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b — > 57 and the measurement of the anomalous magnetic moment of the muon. We 
apply : 

Relic density: 

Recent evidence suggests jSl] that Q x ~ 0.3 with h 2 ~ 0.5. We will take as a 
conservative favored region 

0.03 < fl x h 2 < 0.3. (3.13) 

In addition to our conservative dark matter limit, we also take into account the 
recent results of WMAP 55 that give a 2a range for the density of cold dark matter, 
^cdmIi 2 = 0.1126±°$g. Let us stress that the requirement of (|3.13|) should not be 
treated as an exact constraint, but rather as an indication of the region preferred by 
cosmological considerations because of the uncertainty in mass spectrum calculation. 

b —* 57 Constraint: 

Another observable where the SUSY particle contributions can be important and 
measurable is the flavor changing decay b — > 57 |56j . In the Standard Model, this 
process is mediated by virtual isospin +1/2 quarks and T^-bosons. In supersymmet- 
ric theories, the spectrum allows new contributions involving loops of charginos and 
squarks or top quarks and charged Higgs bosons. For our analysis, we use the results 
given by the CLEO and BELLE collaborations [2ZJ. We adopt the procedure taken in 
the recent benchmark study of Battaglia et al. |58| and choose to impose the constraint 

2.33 x 10~ 4 < BR(6 -> sj) < 4.15 x 10" 4 . (3.14) 
The Muon Anomalous Magnetic Moment: 

Recently, the Brookhaven collaboration has given a new measurement of the anoma- 
lous magnetic moment of the muon |62| 

/ exp _ n\ 

^L- L = a c * p = 11 659 202 (14) (6) x HT 10 , (3.15) 

Following [HE] we introduce the parameter 5^ to quantify the difference between theo- 
retical and experimental determinations of a^: 

Sft = (<V - 11 659 000 x 10~ 10 ) x 10 10 . (3.16) 

From this the current experimental determination of the parameter (L is 5^ xp = 203±8. 
In our discussion, we are less conservative than the authors of jHS] and consider a 2 
standard deviation region about the anomalous moment of the muon based on the r 
decay analysis [5T?] : 

-11.6 < <jnew physics = (5 exp_ J SM < 3Q 4 [2fj] _ (3 ^ 
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4 The Models 



4.1 The general moduli-dominated case 

We present in Figs ^El the parameter space allowed by experimental constraints and 
the corresponding values of two important observables in the context of dark matter : 
the spin-independent scalar cross-section c^ l p and the muon flux from the Sun. 

Let us start with a general comment concerning the relic density. As we can see 
looking back at (|2.1|) and (|2.3j) , the main feature of the moduli-dominated regime is to 
have gaugino and scalar masses mediated by one loop corrections and threshold effects. 
This implies a phenomenology with light neutralinos and possibly light squarks. In any 
case, it has been shown in ^3] that the lightest neutralino x keeps a wino nature in 
a broad region of parameter space, being degenerated with the lightest chargino xt ■ 
The first effect of this degeneracy is a complete depletion of the relic density which 
barely reaches 10 -2 due to the strong coannihilations channel (xxt)- The only way of 
splitting these two masses is by the influence of the universal negative Green-Schwarz 
counterterm Sqs i n <|2.1jl . Increasing \5gs\ decreases the ratio M\jMi leading first 
to the critical value M\\i. e jM2\i. e , = 1 around which the bino and wino contents 
give an interesting relic density. By increasing further I^gsI) we can reach a bino- 
like region, where the relic density is enhanced up to the point where the neutralino 
mass reaches the lightest stau mass giving a density compatible with the last WMAP 
results (f coannihilation corridor). This is very well illustrated on the regions denoted 
O.K O < 0.3 and WMAP in Figs. EJandU 

Figures Hand El (left) present the space constraints in a (t, M 3 / 2 ) plane for Sqs = 0, 
p = and tan (3 = 35. With this choice of parameters, the lightest neutralino is always 
wino-like, independently of the value of t = (T). Indeed, the gaugino masses M a in 
12.11 are proportional to their respective beta-function coefficients b a , which yields the 
standard relation, at the GUT scale : 

M 2 b 2 5 v ; 

implying M 2 < M\ at the electroweak scale. This leaves the lightest neutralino x an d 
chargino %i 111 a wino state, which depletes the relic density of the neutralino (1CP 2 
at its maximum level value). 

If we decrease the value of t from its self-dual point t = 1, the function G2(t,t) will 
reach a point where Mi=i t 2 = 0, corresponding to 

(t + t)G 2 (t,t) = 1 -» Re t = 0.523 (4.19) 

giving a null value for gaugino mass terms, and, as a consequence, for the scalar 
neutralino-proton cross section (proportional to m x ). as shown in Fig. ^ After this 



13 



0.2 



0.6 



0.5 



Fi gure 1: The spin independent scalar cross section in the moduli parameter space for Sqs = 

p = 0, in the (i, M 3 / 2 ) plane, for tan (3 = 35 (left) and (p, tan/3) plane for M 3 / 2 — 40 TeV (right). 
Accelerators and cosmological constraints are given for [i > 0. The labels on the black lines correspond to 
the Logio value of a s ^ l p (pb). For a description of the experimental constraints applied, see Section l3~H 




Tan p= 35. 
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Figure 2: Muon fluxes from the sun in the moduli parameter space for Sqs = P = 0, in the 

(p, M3/2) plane, for tan (3 = 35 (left) and (t, tan/3) plane for M 3 / 2 = 40 TeV (right). Accelerators and 
cosmological constraints are given for fi > 0. The labels on the black lines correspond to the Logio value of 
the flux (/i km" 2 yr" 1 ). For a description of the experimental constraints applied, see Section EOl 
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Tan P = 35. 
t = 0.25 
p = 0. 



Fi gure 3: The spin independent scalar cross section in the moduli parameter space in the 
(Scsi M 3 / 2 ) plane, for t = 0.25, p — 0, /i > 0, for tan (3 = 5 (left) and 35 (right). The labels on the black 
lines correspond to the Logio value of cr*™, (pb). For a description of the experimental constraints applied, 
see Section l3~H 





Figure 4: Muon fluxes from the sun in the moduli parameter space in the (Sgs'i M3/2) plane, 
for t = 0.25, p — 0, /i > 0, for tan (3 = 5 (left) and 35 (right). The labels on the black lines correspond 
to the Logio value of the flux {^1 km~ 2 yr _1 ). For a description of the experimental constraints applied, see 
Section IO 
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pole, Mj = i 2 becomes negative, and its absolute value increases with t : this restores 
a non-vanishing cross-section. This cross-section falls with t becoming small because 
all the sparticles (even the scalars driven by M3 through the renormalization group 
equations) become more massive (G2 — ► 00 as t — > 0), which increases the virtuality of 
the exchanged squarks or neutral Higgs in the elastic diffusion process. Generically, for 
a wino-like neutralino, neutrino flux coming from the Sun are small although annihila- 
tions into gauge bosons (which give more energetic neutrinos than other annihilation 
channels) are favoured. Indeed the capture is small due to suppression of the Z ex- 
change, proportional to the higgsino content of the neutralino. This explains the small 
values of the fluxes in the (t, M 3 / 2 ) plane. This is illustrated on the left panel of Fig. 
121 where muon fluxes (which we will denote by flux®) coming from the Sun are smaller 
than 10 km -2 yr _1 , well below possible detection. 

We show on the right panel of Figs. ^ an d[21the effects of the parameter p and tan (3 
for a fixed value of the gravitino mass (M 3 / 2 = 40 TeV) keeping 5qs = 0. From 1)2.3)) . 
we see that increasing p up to one decreases scalar soft masses, increasing the scalar 
cross section a x °^ p : for p=0.95, the scalars are sufficiently light 1)2.3)1 to allow a cross 
section of the order of 10~ 8 pb in some region of the parameter space. For high values 
of tan (3, niH decreases and one has lower values of \x so higher higgsino content of the 
neutralino before the "noEWSB" boundary which also enhances the Higgs coupling in 
(proportional to 2 x ,i=i,2%,i=3,4) • Along this boundary, although it is excluded by 
limits on b — > S7 and m^, fluxes coming from the Sun for neutrino indirect detection 
can be high. For small tan/3, the cr x c ^ p enhancement comes from the light Higgs 

contribution \q — ► xi- 

Whereas direct detection, being driven by Higgs (h, H) exchange has a strong tan (3 
dependance as can be seen on Fig. ^ the indirect detection {Z exchange in the capture) 
is clearly independent of tan (3 as shown in Fig. [21 

Figures |3] and 0] show, respectively for direct detection and indirect detection, the 
{5qs, M3/2) plane for t = 0.25, p = and tan/3 = 5 or 35. Because of the positive sign 
of 61, the effect of 5gs (which is negative) allows the neutralino to be bino-like and 
splits m x and m x +- Indeed, looking back at 1)2.1)1 . we can see that, at grand unification, 

102 = abl-Scs ' w ^ n a independent of the gauge group considered. Decreasing from zero 
5gs decreases the wino (versus the bino) content of the neutralino and thus increases the 
relic density. For 5qs ~ —23/5, we have the relation Mi\i_ e jM2\i. e . ~ 1 at low energy 
giving £lh 2 \ x ~ VLh 2 \^^j^ p . Taking larger \5gs\ values means driving to the universal 
case, with M\\gut = M?\gut at the high energy scale, and so Mi|^ e ~ O.QM2\i. e . 
at the electroweak scale, giving back a bino-like lightest neutralino. Cosmologically 
favoured relic densities then come from xt coannihilation along the f LSP boundary 
line. However in this region, fluxes are very small as capture is disfavoured by squark 
(though quite light) exchange and very small higgsino content of neutralino in Higgs 
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M3/2=5 TeV 




Figure 5: The spin independent scalar cross section as a function of the Green-Schwarz coun- 
terterm Sqs, for t = 0.25, p = and different values of M 3 / 2 , with tan j3 = 5 (left) and 35 (right). 
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Figure 6: Muon fluxes from the sun as a function of the Green-Schwarz counterterm Sgs> for 
t = 0.25, p = and different values of M 3 / 2 , with tan j3 = 5 (left) and 35 (right). 
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Moduli case vs Direct Detection : no exp cut 



Moduli dominated SUSY breaking vs Direct Detection 




Iv^ (GeV) M z (GeV) 



Figure 7: The spin independent scalar cross section as a function of the neutralino mass M x 

for t — 0.25, p — and tan/3 = 5, 20, 35. We have scanned the moduli parameter space on M 3 / 2 and Sqs 
before (left) and after (right) having applied the accelerator and cosmological constraints of Section l3~^l 



exchange. 

Decreasing 5gs first decreases IM3I down to gluino LSP points and then increases 
IM3I (which is now negative), scalar masses (through M3 RGE effect) and \i (decreasing 
the neutralino higgsino fraction). In particular the lightest Higgs h (through its stop 
radiative correction) and the heavy neutral Higgs H masses increase and the coupling 
z X:i(2) z x.3(4) decreases so that a x c ® l p decreases following the higgs mass experimental 
limit contour (see Fig. |3J). For high values of M 3/ / 2 , ^x-p * s ac t uau y beyond reach 
of experiment sensitivities (< 10~ n pb). Moreover, for small values of SqSj a x-p 
follows the contour of iso-^ values (the higgsino content of neutralino) along the m + 
limit. Higher values of tan/3 are also more favourable to direct detection because the 
Higgs H is lighter. For low values of M3/2, \i values are smaller and the higgsino 
fraction of neutralino drives the phenomenology. Relic density is closed to the WMAP 
range and both direct [p^L ~ 1CP 8 pb in Fig. EJ) and neutrino indirect detection 
(flux® ~ 10 3 ~ 2 km _2 yr _1 in Fig. can be interesting but this region is excluded by 
limits on and b — * 57. Muon fluxes coming from the Sun follow the iso— fi shape 
given by the limit on m +. Small \Sgs\ can a ls° l ea d to gluino LSP. 

We have illustrated explicitly the 5gs dependance of a x c ^ l p and flux® on Figs. El 
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Moduli case vs neutrino Indirect Detection : no exp cut 



Moduli dominated SUSY breaking vs neutrino Indirect Detection 




700 
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Figure 8: Muon fluxes from the sun as a function of the neutralino mass M x for t = 0.25, p = 
and tan/3 = 5, 20, 35. We have scanned the moduli parameter space on M 3 / 2 and Sqs before (left) and 
after (right) having applied the accelerator and cosmological constraints of Section ITU 



and|21 As previously said, decreasing 5gs increases scalar masses (through M3 running 
effect) so m,H which dominantly drive a^ p . It also exchanges wino and bino content 
leading to smaller coupling (taxiOw suppressed) and lower values of a^_ p . We also 
see the well known strong influence of tan/3 on 771/7(771^) leading to higher values for 
tan/3 = 35. 

Concerning indirect detection, bino content decreases the fluxes but small values 
°f \$gs\ can lead to gluino LSP and high fluxes nearby. This corresponds to the effect 
of small values of M^\qut in soft parameter running, leading to small values of \x and 
favoring both direct and indirect detection [S]. But in the models we consider in this 
section, this happens in region excluded by accelerator constraints because of light 
chargino contribution in b — > S7 and too light Higgs h. 

To conclude the discussion of the moduli-dominated SUSY breaking case, we show 
a wide range of models in the (m x ,a^: p ) plane for direct detection (Fig. 0) and in 
the (m x ,flux®) plane for neutrino indirect detection (Fig. |HJ). We then clearly see 
the complementarity between dark matter and accelerator searches by comparing left 
panel (no accelerator constraints) and right panel (accelerator constraints applied) of 
these two figures. 
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In regions satisfying accelerator constraints, moduli dominated model give rise to 
either wino neutralino with strongly suppressed relic density and very low detection 
rates, or bino neutralino (thanks to Sqs) an d interesting relic density via xt coanni- 
hilation or adjusted wino vs bino content but again the direct and indirect detection 
rates are beyond reach of detection (except sometimes for a larger than one ton size 
direct detection experiment such as Zeplin Max). This situation is summarized on 
Figs. El and El 

One interesting feature of this scenario is that Mi—\ t 2 is becoming negative for 
large values of |5gs|- This negative sign can allow cancellation between the up-quark 
and down-quark contribution to the scalar cross section o*£- p together as for negative 
values of ji pointed in jSU]. This phenomenon is responsible of the "seagull" shape of 
the plots El and depleting in the same way the scalar cross section down to 1CT 13 
(for tan/3=5) or 10 -11 (for tan/?=35). 

Moduli models satisfying accelerator constraints and a WMAP favoured relic den- 
sity are generically not detectable by dark matter searches. 

4.2 Dilaton dominated case 

The phenomenology of the dilaton dominated scenario is completely different from 
the moduli domination just discussed. If we look at Eqs <|2.7|) and 1)2. 5 J) , it is clear 
that we are in a domain of heavy squarks and sleptons (of the order of the gravitino 
scale) and light gaugino masses (determined by the dilaton auxiliary component vev). 
Indeed, the beta-functions b a being of the order of 10 -2 , the corresponding terms are 
not competitive by comparison to the F term of the dilaton in 1)2. 5 Jl . In fact, looking 
more closely at 1)2.9)1 and (|2.8jl . for not too large values of the universal beta-function 
coefficient of the first condensing group (&+), we may consider that F s is a linear 
function of 6+. Increasing 6 + means approaching the universal case for the gaugino 
sector (and the scalar one, driven by M 3 / 2 ). 

Figures l9l andlTOl present the (b+, M 3 / 2 ) plane for tan (3 = 5 and 35 with experimental 
exclusions, neutralino relic density and respectively iso-cr®^ and iso-flux® curves. We 
first discuss the experimental exclusion plots. 

For low values of M 3/ / 2 , scalars and gauginos are light so that accelerator constraints 
are strong. At fixed values of M 3 / 2 and decreasing values of b + , we see from (|2.7|) that 
gaugino masses decrease 2 . Ms\gut being smaller, the m 2 H running slope is softer and 
yields positive mfy at low energy leading to lower values of This explains the region 
with too light a chargino (i.e. Higgsino) mass followed by the "no EWSB" region as 
one goes along a decreasing b + direction. 

2 See Fig. For very low values of 6+ around 0.1 (independent on M 3 / 2 value), we have a very localised 
range with gluino LSP due to cancellation in M3 coming from the cancellation between 63 (< 0) and 6+. 
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As one increases M 3 / 2 at fixed b+, one decreases M 3 and the same focussing effect 
as discussed in the previous paragraph leads to a lower value of the SUSY parameter 
fj,. It is thus not surprising that the "no EWSB" extends further to the right as M 3 / 2 
increases. 

We now turn to a discussion of the direct and indirect (neutrino) rates. For low 
values of the gravitino mass, these rates are favoured because of the light squark con- 
tributions in v^-p and in 0"5^ (enhancing capture and flux®). 3 

As one increases M 3 / 2 , one finds a region where the lightest neutralino is a mixed 
higgsino-gaugino state. It can satisfy WMAP requirements on relic density through 
XX ~^ W + W~ , ZZ or ti annihilation and XXi >XX% coannihilation processes. This 
region is safe from limits on rrih and b —* 57 thanks to the high M 3 / 2 values. Direct 

detection is favoured because the higgsino component enhances coupling in XQ — * XQ 
leading to v^-p ~ 10~ 7 ~ 8 pb. In this mixed region, indirect detection is also powerful 

(Fig. I1U|) . Indeed, the higgsino component allows efficient capture via XQ — * XQ an d 
neutralino annihilations into gauge bosons or ti lead to energetic neutrinos/muons. 
Muon fluxes coming from the Sun are enhanced and can reach 10 1 to 10 3 km~ 2 yr _1 . 
In both cases (direct or indirect detection), the detection rates increase with tan/3 
because the higgsino fraction is higher, direct detection being also enhanced by a 
lighter H Higgs. 

The interesting effect of b + on both direct (er*™,) and indirect (flux®) detection 
is illustrated on Fig. 1121 We see that decreasing 6 + can give a gain of one order of 
magnitude for direct detection and up to 2 orders of magnitude for indirect detection 
which is favoured twice by the higgsino fraction (capture plus annihilation in gauge 
bosons giving more energetic neutrinos/muons). cf^- p being proportional to z Xj iz Xj3 ( 4 ), 
finally decreases when 6+(/i) is further decreased (% too much higgsino like). The bumps 
on Fig. El (left panel) occur around the [i-M\ crossing (Fig. Ill|) . This requires higher 
values of b + for higher values of M 3 / 2 . The last re-increasing of cr^ l p (for M 3 / 2 = 1 TeV) 
comes from the \i evolution (Fig. HT|) and Zx> lZ xM4) dependence. 

Fluxes® have the same behaviour with M 3/ / 2 and &+, but being proportional to 
z x3(4)' wnen ^+ decreases they still increase or actually follow the (inverse) evolution 
of the parameter /i as can be clearly seen by comparing Figs. ^2 and El (right panel). 
Starting from high values of b + , we can also notice the first bump corresponding to 
m x = mt because neutrino coming from XX ~> annihilation are less energetic than 
XX ~> W + W~ , ZZ. This very favourable effect on neutralino dark matter relic density 
and detection rates consisting in decreasing \i and scalar masses through M^gut run- 
ning effect has been pointed out in |5j but is directly related here to the 6+ parameter. 

To conclude this dilaton case, if we compare model predictions with experimental 

3 b + acts in the same way since we have seen that it is increasing with /j,. But these low M 3 / 2 regions are 
ruled out by accelerator constraints. 
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Tan p=35 



Figure 9: The spin independent scalar cross section in the dilaton parameter space, in the (&+, 
M3/2) plane, for tan (3 — 5 (left) and tan (3 = 35 (right). Accelerators and cosmological constraints are given 
for /i > 0. The labels in the black lines correspond to the Logio value of c£™p (pb). For a description of the 
experimental constraints applied, see Section 




Figure 10: Muon fluxes from the sun in the dilaton parameter space in the M 3 / 2 ) plane for 
tan/3 = 5 (left) and 35 (right). Accelerators and cosmological constraints are given for /1 > 0. The labels in 
the black lines correspond to the LogiQ value of the flux (p km~ 2 yr -1 ). For a description of the experimental 
constraints applied, see Section EOl 
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Figure 11: Low energy parameter M 1 , M 2 , M 3 , fi evolution with b + , for tan/3 = 35 and M 3/2 =1000, 3000 
GeV. 



flux© ~ 



sensitivities for both direct (Fig. left panel) and neutrino indirect detection (Fig. 
1131 right panel), we see that dilaton models satisfying both accelerator constraints 
and approximate WMAP relic density give mixed higgsino-gaugino neutralino with 
generically high detection rates : tt x h 2 ~ (nh 2 )^lf p , a s x ™ l p ~ 1(H 7 to 9 ) pb and 

to 3]^ m -2 yi--l_ 

One of the most interesting feature of the dilaton type models is its closed parameter 
space. Indeed, the beta function of the first condensing gauge group cannot exceed the 
largest one of the models ipES ~ 0.57). On the other way, the No EWSB condition 
forbids high values of M 3 / 2 , the upper limit depending on tan/3. Imposing WMAP 
constraints on the relic density, in this this closed (M 3/ / 2 , 6+) plane, neutralino mass is 
limited to 500 GeV (resp. 1500 GeV) for tan/? = 35 (resp. 5). Those limits can be 
compared to the mSUGRA case where m x < 500 GeV for tan/3 < 45 3 (small values 
of Mq and M]y 2 ) but where this bound can be higher (~ TeV) if one considers the 
hyperbolic branch/focus point region 4 . 

Dilaton models could be detected by future dark matter searches especially a ton- 
size direct detection experiment like Zeplin [SJ, or, for points with m x < 500 GeV a 
km 3 neutrino telescope like Icecube 
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Figure 12: The spin independent scalar cross section (left) and the muon flux from the sun 
(right) in the dilaton parameter space, as a function of b+ for tan /3 = 35 and different values of M 3 / 2 



5 Conclusion 

In the specific context of the class of string models that we have considered, we have 
seen that the predictions regarding dark matter are strikingly different according to 
the type of supersymmetry breaking considered. In the case of moduli domination, one 
does not expect any signal in the forthcoming direct or indirect (neutrino) detection 
experiments. On the other hand, these experiments should not miss the neutralino 
signal in the case of dilaton domination. Thus the detection of dark matter or the ab- 
sence of detection may give key information on the nature of supersymmetry breaking, 
at least in the context of this given class of models. 

Obviously there are connections between these results and detection of the LSP at 
colliders. Small direct detection cross sections or small indirect detection fluxes are 
obviously correlated with small production cross sections at colliders. In any case, 
it is interesting for collider searches to note the characteristics of the regions that 
satisfy the criterion of satisfactory relic density. For moduli domination, we have 
identified two regions of interest: one where m v ~ m + ~ m^o through the bino 
and wino content of the LSP (for sufficiently large values of tan/3), and the other one 
close to stau LSP region where m^ x ~ m x . In the case of dilaton domination, the 
cosmologically interesting region corresponds to a LSP with a proper higgsino content 
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Dilaton dominated SUSY breaking vs Direct Detection 



Dilaton dominated SUSY breaking vs neutrino Indirect Detection 




Iv^ (GeV) M z (GeV) 



Figure 13: The spin independent scalar cross section (left) and the muon flux from the sun 
(right) in the dilaton parameter space, as a function of the neutralino mass M x after a scan on b + 
and M3/2 for tan j3 = 5, 20 and 35. Accelerators and cosmological constraints (taken for /1 > 0) have been 
included. For a description of the experimental constraints applied, see Section 



{m x ~ m x + ~ m x°)' Furthermore the parameter space being closed, this case gives an 
upper bound on neutralino mass: m x < 1500 GeV. 

We think that the type of conclusions that we have reached is relevant to a larger 
class of models. It remains however to make similar analyses on other classes of models 
to see if dark matter searches, whether they are positive or negative, give interesting 
indications on the way supersymmetry is broken, and on the way it is transmitted to 
the observable sector. 
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